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[ — ■ Abstract We discuss the e-method as used in various recent QED 

bound-state calculations by considering mathematical model examples. 
Recently obtained results for higher-order self-energy binding correc- 
tions at the two- loop level are reviewed. Problems associated with the 
interpretation of squared decay rates as radiative bound-state energy 
level shifts are discussed. We briefly expand on the relation of squared 
decay rates to nonresonant and radiative corrections to the Lorentzian 
_C ■ line shape, including their dependence on the experimental process un- 

der study. 

X 

.«?. 

PACS numbers: 31.15.-p, 12.20.Ds 

Keywords: Calculations and mathematical techniques in atomic and molecular physics, 
quantum electrodynamics - specific calculations. 



CO 

o 



1 



1 Introduction 



This paper is concerned with mathematical methods employed in recent analytic evaluations 

El El °f higher-order binding corrections to the Lamb shift. These methods rely on a 
separation of the virtual photon energy integration into high- and low-energy domains. The 
methods are applicable in a wider context, and we attempt to provide a certain clarification 
by considering mathematical model examples. We focus on the two- loop self-energy correction 
(see Fig. ^) in hydrogenlike systems with a low nuclear charge number. We also discuss related 
corrections in helium. 

In the second part of the paper, we discuss a recent investigation [2J which is concerned with 
predictive limits of energy shifts as derived from the Gell-Mann-Low-Sucher Theorem |S1E]. 
Expressions obtained based on this theorem have provided the basis of level-shift calculations for 
decades; these may not be accurate enough for projected future experiments. Certain problems 
associated with this theorem find a rather natural solution in the two-time Green function 
method |1U| . other problematic aspects connected with this theorem concern the interpretation 
of level shifts involving squared decay rates. 

Our investigations are motivated by the recent dramatic progress in laser-spectroscopic experi- 
ments in atomic hydrogen (e.g. [HJESEIl) as well as helium [miTClUHllTTj . 



(a) (b) (c) 

Figure 1: The crossed (a), rainbow (b) and the loop-after-loop diagram (c) which 
contribute to the two-loop self-energy for a bound electron. The propagator of the 
bound electron is denoted by a double line. 



2 The "Epsilon Method" 

In QED bound-state calculations, we are often faced with the following problem: how to separate 
terms of a given order in the (Za)-expansion, and how to treat infrared divergences that originate 
from higher-order terms in the expansion of the bound-electron propagator in powers of the 
binding field. 

The so-called "e-method" has been employed in the analytic calculation of self-energy effects 
in bound systems. The energy scales for the virtual photons are treated separately: (i) the 
nonrelativistic domain, in which the virtual photon assumes values of the order of the atomic 
binding energy, and (ii) the relativistic domain, in which the virtual photon assumes values of 
the order of the electron rest mass. The two energy domains are separated via a parameter 
e. Without appropriate approximations and expansions applicable to the two energy domains, 
respectively, the analytic evaluation of either the high- or the low-energy part would not be 
feasible. 

In one-photon calculations, we have to deal with one virtual photon energy u. For two-photon 
problems mElE], one has to generalize the method to the case of two virtual quanta and, by 
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consequence, two separation parameters t\ and €2 (see also Fig. 1 of In both the one- and 
the two-photon case, we require the dependence on the separation parameters to vanish at the 
end of the calculation, i.e. after the high- and the low-energy parts are added. 

We follow here the discussion in App. A of Ref. jS], and we consider a model problem with 
only one "virtual photon". In contrast to Ref. [Sj, we choose a mathematical model problem 
of a slightly more complex structure in order to illustrate the occurrence of double-logarithmic 
terms in the semi- analytic expansion, which involves powers and logarithms of the expansion 
parameter. The mathematical model example reads 

J(P)=j\n(co) ^t^du,, (1) 

where the integration variable u might be interpreted as the "energy" of a "virtual photon". 

We intend to derive a semi-analytic expansion of J(f3) in powers of [5 and In /3. The quantity Za, 
which parameterizes the strength of the binding Coulomb field, replaces the expansion parameter 
(3 in actual self-energy calculations. The "high-energy part" of the calculation is given by the 
expression 

JnW) = J\(lj)^^6lj. (2) 
In the integration domain u> £ (e, 1), we may expand 



V^+^ = u; + ^ + ^ + 0(/3 6 ). (3) 
zuj our 

However, this expansion is not applicable in higher orders to the domain to £ (0, e) because of 
the appearance of inverse powers of lv which lead to "infrared divergences" . 

We expand the integrand of Jh(/3) first in powers of (3 according to Eq. ©. The resulting integrals 
can be evaluated analytically. Every term in the /3-expansion is then expanded in powers of e up 
to the order e°. Higher-order terms in e are irrelevant; they cancel at the end of the calculation 
(just as the divergent terms in e), because the original expression for the integral J(/3) in Eq. Q 
is manifestly independent of e. The result of the calculation of Jn(f3) is 

MM = {ln(2)-l + 0(e 2 lne)}+/? 2 |-I ln 2 ( e ) + ^ ln 2 (2) - ^ + 0( e 2 lne)} 

+(3 4 I— In 2 (e) - — ^ In (e) - - — In 2 (2) 
\ 32 W 16 e 2 w 32e 2 32 v ' 

+— In (2) + — + — + 0(e 2 In e) 1 + 0{(3 6 In 2 (3) . (4) 
32 v ; 384 64 v 1 J KH H> y J 

The contribution Jh(/3) corresponds to the "high-energy part" in analytic self-energy calcula- 
tions, where the propagator of the bound electron may be expanded in powers of Za [see Sec. Ill 
of@]. 

The expression for the low-energy part uj £ (0, e) reads 

J h ((3) = I ln(w)y T f£-da;. (5) 



3 



We have to keep the numerator of the integrand \J lo 2 + (3 2 in unexpanded form. However, within 
the integration domain uj G (0, e), we may expand the denominator \/I — uj 2 of the integrand in 
powers of uj. Because uj < e, the expansion in powers of u is in fact an expansion in (3 in the 
low-energy part. 

One may draw an analogy between the term \/uj 2 + f3 2 and the Schrddinger-Coulomb propagator 
in self-energy calculations m^HSHS]- In the low-energy domain, this propagator may not be 
expanded in powers of the binding field. The expansion 

_=^ 1 + T + ^ 1 + ^> (6) 

corresponds to the (Za)-expansion in the low-energy part. In actual self-energy calculations (a 
detailed discussion can be found in Ref. [Hj), higher-order terms in the low-energy part originate 
from the generalized Foldy-Wouthuysen transformation of the transition current, from relativis- 
tic corrections to the Foldy-Wouthuysen transformed Hamiltonian, and higher-order terms in 
the multipole expansion. Specifically, more details concerning the multipole expansion can be 
found in the discussion following Eq. (11) of Ref. 

We obtain for the low-energy part, 

MM - ✓{itfW-i h,(|) -i-^ + ofm.)} 

+ ^{-s h> w + ii? b w + s? + s h, (!) 

--ln(^) + — - — + 0(e 2 lne)l +0(/3 6 ln 2 /3). (7) 
64 \(3 J 192 256 v ') v 1 y J 

When the high-energy part (jlj) and the low-energy part (JJJ) are added, the dependence on e 
cancels, and we obtain the result 

J {13) = J H (/3,e) + J L (/5,e) 

= ln(2) - 1 

This result clearly demonstrates the semi-analytic character of the expansion: it involves double- 
logarithmic terms In 2 (2/(3) and single logarithms In (2/0) as well as constant terms. The same 
pattern is observed in actual self-energy calculations. 
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3 One— and Two— Loop Self— Energy Calculations 



It is well known that the one-photon self-energy correction to the binding energy in \ow-Z 
hydrogenlike systems can be parameterized as as 

where the dimensionless quantity F(Za) has the following semi-analytic expansion. 

F{Za) = An ln(Za)- 2 + A 40 + (Za) A50 

+(Za) 2 [A32 ln 2 {Za)- 2 + Asi ln(Za)" 2 + Ago + K] , (10) 

where 7Z vanishes as Za — * 0. The Acoefficients are state-dependent. In the following, we focus 
on P states. The coefficients A41 and A52 vanish for P states and states with higher orbital 
angular momenta. 

In order to illustrate the analogy of our mathematical model example (see Sec. |2J) with one- 
loop self-energy calculations, we give here the high-and low-energy parts derived in [2] for the 
self-energy of an electron bound in a hydrogenlike system (the 2Pi/ 2 state): 



1 ,2 



F}1 (2P 1/2 ) = -- + (Za) 



6 

and 



2 103, , . 4177 103, , n . . 

In e + In 2) + 0(e) 

9e 180 w 21600 180 w w 



+ 0(Zaf, (11) 



F L (2P 1/2 ) = -^lnA ; o(2P) + (Za) 2 



l + Wo l ^TzW ] -°- 79mi)+0 ^ 



+ Q{Zaf . 



(12) 

Adding the two contributions, the dependence on e cancels, just as we had observed when forming 
the sum of the high-energy part and the low-energy result Q. 

The two-loop self-energy correction to the energy of hydrogenlike systems with low nuclear 
charge number reads 

6E SE = (Zar^H(Za), (13) 

where the dimensionless function H(Za) is given by 
H{Za) = B 40 

+(Za) 2 [B 63 ln 3 {Za)- 2 + B m \n 2 (Za)- 2 + B 61 ln(Za)~ 2 + B m + W] , (14) 

where 1Z' vanishes as Za — > 0. In two-photon calculations, we introduce two separation pa- 
rameters ei and This leads to four different integration regions: (i) both photon energies 
are small (the "low-and- low-energy part"), (ii)+(iii) two mixed contributions (with one large 
and one small photon energy), and (iv) a "high-and-high-energy part" with two large photon 
energies. 

Bq3 vanishes for all P, D, F, G, . . . states, i.e. for all atomic states with a nonvanishing orbital 
angular momentum, and we first discuss here the the coefficient Bq 2 - The low-and-low-energy 
part of the two-loop problem (both virtual photons have a small energy) can be evaluated 
using nonrelativistic quantum electrodynamics (NRQED). The relevant expression is given in 
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Eq. (16) of Ref. The following double-logarithmic term of order a 2 {Zaf 1 originates from the 
low-energy part [see Eq. (38) of Ref. ..6\], 



In 



ei 



(Za) 2 



In 



f2 



{Zaf 



2Tr(A5^(r)) 



9 m 4 



where the known result for the matrix element reads 
(A5 i3) (r)) = A 



'Jn,l=l,m \ 



r=0 



(15) 



(16) 



Because Bq% vanishes for P states, this result for the low- low-energy part determines uniquely 
the total result for B§2- This is because the dependence on e\ and e 2 necessarily has to cancel 
at the end of the calculation according to 



In 



{Za) 2 m 
+ ln 



In 



02 



In 



(Za) 2 m 
(2 



, 1 m 
+ ln ( — 

£2 



In 



(Za) 2 m 



+ In 



111 



ei 



In 



(Za) 2 m 
m 

(2 



ln 2 [(Za)" 



(17) 



1,2), 



Note that the logarithm ln[ei/(Za) 2 m] is characteristic of the low-energy domain 
whereas the logarithm ln(m/ej) is characteristic of infrared divergencies that occur in the eval- 
uation of integrals involving highly energetic virtual photons. In view of Q17|). we may conclude 
that the coefficient of 



In 



ei 



(Za) 



In 



£2 



[Zaf 



in the low-and-low-energy part is the same as the total coefficient of ln 2 [(Za)~ 2 ] for the two-loop 
self energy. This has lead to a rigorous derivation of B§2 for P states 6 , confirming the results 
of the previous investigation ^ 



4 n 2 — 

5 62 (nP) = 5= — 
27 n z 



1 



(18) 



This result is valid for all P states independent of the electron spin. 

The e-method provides a convenient tool for the analysis of the problematic nonlogarithmic 
corrections an d Bqq. Let us recall that the evaluation of the one- loop coefficient f° r S 
states has attracted attention over many years |19U2Ull2~HI22l lT]. Today, we can hope to evaluate 
the corresponding i?60-coefficients using this method. 

As a first step in this direction, we have obtained results [3] for the following fine-structure 
differences of -E>6fc-coefficients (fc = 0, 1) of P states, 



A is B 61 = B 6 i(nP 3 / 2 ) - B 61 (riPy 2 ) , 
A fsJ B 60 = B 60 (nP 3/2 ) - B m (nP 1 / 2 ) . 



(19) 



We implicitly define the symbol Af s to denote the difference of the coefficient for an nP 3/ / 2 state 
and the corresponding nPi/2 level. Note that the fine-structure difference of the leading double 
logarithm vanishes [see Eq. (|18|)]. 

A fsJ B 62 = . (20) 

Certain simplifications are possible when considering the fine-structure difference of the Bqi and 
-B60- coe ffi c i en ts. Specifically, the contribution of the high-and-high-energy integration domain 
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can be investigated with the help of a modified Dirac Hamiltonian. In this context, vertex 
corrections are taken into account by considering the electron form factor F\ and F 2 . A further 
simplification occurs because it is possible to devise a unified treatment for both the low-and- 
low-energy domain and the mixed-energy contributions. Some of the mixed-energy effects can 
be described by magnetic form-factor corrections to the leading one-photon self-energy. Because 
the magnetic form factor F 2 does not have infrared divergences, contributions of the type 



In I - I In 

e, 



£3- 



(Za) 2 m 



(i = l,2) 



vanish for the fine-structure difference. We are therefore left with an infrared-safe and (in the 
context of the effective form-factor treatment) also ultraviolet-safe mixed-energy integration 
domain, for which a simplified treatment is possible. 

We recall here the relevant results from [3] . The high-and- high-energy integration domain yields 



Er = Ei + E 2a + E 2 b + E 3 , 



(21) 



where the correction E\ is due to the F\ form factor, E 2a and E 2 b are due to the magnetic F 2 
form factor, and E3 is caused by a second-order effect involving two one-loop magnetic form 
factor corrections to the spin-orbit interaction. The results read Eq. (23)] 



Ei 



a\ 2 (Za) 

71 



6 



n° 



F{ {4) (0) 



n 



1 



73, z 



773 



(22) 



and according to jSJ Eq. (27)] 



\7r/ n 6 



F^rm ( m ,& 23 \ 

F2 (0) \T2n + ^-m?) 



m , 



(23) 



as well as Eq. (36)] 



E 2 b 



a\2 {Zaf 
ttJ n 3 



1 n 2 



6 n 2 



m m 
In h In 

2ei 2e 2 



+ 2^ (4) (0) m 



18 



n 



73. 



773 , (24) 



and according to |SJ Eq. (30)] 



£ 3 =(- 

v 7T 



{Zaf 



TV 



227 1 
+ 



2880 3273 80n 2 



773 . 



(25) 



Analytic results are known {2312412 5 for the two-loop form factors entering into these expressions 
(we take into account only the diagrams in Fig. ^ and ignore the vacuum-polarization insertion 
in the one- loop vertex correction), 



T^F^O) 



•— -— C(2) + 3C(2) In2--C(3), 
576 144 SV ; SV ' 4 SV ; ' 



(26) 



F?\0) 



773 2 ^ (4) (0) 



315 3 
.£i + | C (2)-3C(2) ln2 + |c(3), 



•^ + ^C(2)--C(2) In2 + ^C(3) 
240 40 sw 10 sw 40 sv ' 



(27) 
(28) 



7 



The sum of the low-and-low-energy domain and the mixed integration regions is 

Eh = E 4 + £5 , 
where the contribution E^ reads [5J Eq. (51)] 



E 



E^a + E. 



4b — — 



(Za) 6 m 



1 



In- 



2ei 



+ In 



2e 2 



n u [ 6n 2 \ [Zol] 2 m 

and the explicit results for the ^4(72) are given by |3ll26| 

A fs £ 4 (2) = 0.512 559 769(1), 

A fe 4(3) = 0.513 111 333(1) , 

A fs 4(4) = 0.516 095 539(1) , 

A fs £ 4 (5) = 0.518 940 860(1). 



(Za) 2 772 



+ 



22 



1 



22^ 



(29) 

A fs £ 4 (n) , 
(30) 



E 5 reads Eq. (56)] 



where 



E, 



/a\2 (Zaf'm 
\irJ n 3 



72 



1 



A fs 4(n) 



(31) 
(32) 



A fs 4(2) = -0.173 344 868(1) , 

A fe 4(3) = -0.164 776 514(1) , 

A fs 4(4) = -0.162 263 216(1) , 

A fs 4(5) = -0.161 165 602(1) . 



(33) 



Adding all contributions E% 



E$, the dependence on both ei and e 2 cancels, and we obtain 
72 2 - 1 



fs-D61 



3ri 2 



(34) 



as well as 



fs-°60 



( 227 1 
2880 + 32^1 



8O72 2 



+- 



72 



1 



22 z 



+ E 



(4),S, 



487 



(0) — + — 
1 720 877 



23 \ 
2077 2 J 



F{ (4) ' S (0) + 2^ (4) ' S (0) ) 772^ - - + A is h(n) + A fs 4(n 



The explicit results for the principal quantum numbers n = 2-5 read 

A fs 5 60 (2) = -0.361 196 470(1) , 



Af s -B6o(3) 
A fsJ B 60 (4) 
A feJ B 60 (5) 



-0.410 149 385(1) , 
-0.419 926 624(1) , 
-0.420 872 513(1). 



(35) 



(36) 
(37) 
(38) 
(39) 



These results have recently been generalized to the case of helium [27] and used for an estimate 
of higher-order binding corrections to the large and small fine-structure intervals of the triplet 
P levels. 

At this point we would like to mention the further recent progress in the understanding of higher- 
order binding two-loop self-energy corrections (see Refs. 



, which is faced with a number 
of conceptual and calculational difficulties, both in analytic as well as numerical approaches. 
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4 Squared Decay Rates 



An intriguing problem of bound-state quantum electrodynamics is the interpretation of squared 
decay rates which follow from the Gell-Mann-Low-Sucher theorem when applied to excited 
atomic states in two- loop order. We have recently shown (Ref. |7j) that the squared decay rates 
cannot be interpreted in a natural way as real energy shifts. Roughly speaking, the problems in 
the interpretation originate from the fact that the Gell-Mann-Low-Sucher formalism involves 
a priori asymptotic states with an infinite lifetime (vanishing decay rate). The decay rate which 
enters in one-loop order adds an imaginary part to the energy whose square cannot be interpreted 
consistently as an energy shift within a formalism whose starting point was a theory that involves 
asymptotic states with zero decay width (for a more detailed discussion see Ref. [7j) Rather, a 
part of the problematic corrections can be incorporated in a natural way in a modified bound- 
electron Green function according to Eq. (27) of Ref. [7j which involves a "decay-rate operator" 
f 1 defined in Eq. (24) of Ref. 0. Our formula (27) of Ref. 0, which is equivalent to Eq. (16) on 
p. 218 of Ref. |31| . could be interpreted to suggest that the modified Green function simply has 
a pole on the second (unphysical) sheet of the Riemann surface (see also the discussion on p. 217 
of Ref. |B5)- However, this is not the case: the only correction of the "squared-decay rate" type 
which can be incorporated in a natural way into the electron Green function formalism is the 
one caused by the loop-after-loop diagram in Fig. ^ (c), which is discussed in Eqs. (9) - (14) 
of Ref. [7j. The interpretation of the other problematic corrections of the type of a "squared- 
decay rate" discussed in Ref. [see Eqs. (5), (16), and (18) ibid.} cannot be given as easily. In 
order to go beyond the predictive limit of the current theory set by the squared decay, one has 
to consider, in a fully gauge- invariant treatment, the excitation of the atom from the ground 
state via the absorption of (laser) photons, and the return to the ground state via spontaneous 
emission. Resummations of sets of diagrams near resonance may be required. It has been stressed 
in Ref. [3| that the ground state is the only "true" asymptotic state which may be used as an 
in- and out-state in scattering theory. 

We would like to stress here that the two-time Green function method |1U| avoids a number of 
problems associated with the Gell-Mann-Low-Sucher theorem. Aside from the simplified treat- 
ment of degenerate states, we would like to mention the well-known fact that the infinitesimally 
damped S^A-matrix [see Eq. (2) of Ref. [Zj] is, strictly speaking, not renormalizable because the 
damping parameter e breaks the covariance. 

A further problematic aspect of spectroscopic measurements is given by the nonresonant cor- 
rections 32, 33, 3H EE]- It has been stressed in [S3] that nonresonant terms are enhanced in 
differential vs. total cross sections, and estimates for the effect in hydrogenic S-P transitions 
have been obtained (see Sec. 3 of 34 J. The enhancement of nonresonant effects in differential as 
compared to total cross sections also follows in a natural way from the two-time Green function 
formalism [see Eqs. (3.1.120) and (3.1.121) of [Ej]. We observe that the experimental accuracy 
is approaching the 1 MHz level at which the nonresonant terms become relevant (e.g. |36|). 

The order-of-magnitude at which the nonresonant terms enter in two-photon transitions depends 
crucially on the process under study |34U35| . After the two-photon absorption, the atom may 
return to the ground state via spontaneous emission of two photons. In this case, the contribution 
of nonresonant terms (on the level of 10 -14 Hz) is negligible at current and projected levels of 
experiment accuracy |33]- However, this does not imply that the experimental line shape should 
remain Lorentzian up to this level of accuracy. For this process, the dominant correction to 
the Lorentzian line-shape in the two-photon transition is given by radiative (not off-resonant!) 
corrections, and an estimate of a relative contribution of order a (Za) 2 has been given in 
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(see Fig. 3 ibid., this translates into ~ 10~ 6 Hz for the 1S-2S transition in atomic hydrogen). In 
the current experiment JS]; the excited hydrogen atom (2S) is quenched to the rapidly decaying 
2P state. This leads to an experimental line width of the order of a kHz. In this case, the 
nonresonant terms are enhanced, as argued in [221, an d an estimate of nonresonant corrections 
of the order of 10 -2 Hz has been given for this different experimental setup. 

In all cases where off-resonant effects were considered, a formula of the general structure 

[experimental decay width T] 2 



[typical atomic energy level difference AE] 



(40) 



has been obtained for the magnitude of the problematic shift of the peak of the photon scattering 
cross section (which is a nonresonant correction to the Lorentzian line shape). The meaning of 
the "typical atomic energy level difference" depends on the process under study: for differential 
cross sections in hydrogenic S-P transitions, a fine-structure level difference should be used for 
AE (see Ref. |34p. whereas for total cross sections, the correct estimate is obtained by inserting 
energy differences between states with a different principal quantum number (see Ref. [SSI)- 
The order-of-magnitude estimate (|40|) implies that (i) nonresonant effects are smaller than the 
experimental line width by roughly a factor of T/AE and (ii) the magnitude of the nonresonant 
terms decreases at decreasing experimental line width. A formula similar to (|40j) can be used in 
order to estimate the order-of-magnitude of the energy level "shift" by squared decay rates. 



5 Conclusions 

We have discussed the evaluation of higher-order binding corrections to the one-and two-loop 
self energy via the e-method (Sec. |2J) This method has proven to be a useful calculational tool, 
as it leads to a rather clear separation of terms which contribute at different orders in the 
(Za)-expansion and to a transparent formulation of the physical problem. A first step in a 
systematic investigation of the highly problematic nonlogarithmic i?60- coer ficient for hydrogenic 
bound states is presented. The cancellation of the expansion parameter e at the end of the 
calculation is demonstrated by way of a mathematical model example [Eqs. Q - ©] and in 
concrete QED bound-state calculations at the one- and two- loop level (see Sec. EJ). 

When performing two-loop self-energy calculations for excited states, one is lead in a natural way 
to the problem of the interpretation of squared decay rates. These effects cannot be interpreted 
self-consistently as radiative energy shifts of a specific atomic energy level (Sec.@J). We consider 
the connection to the Gell-Mann-Low-Sucher theorem and to nonresonant corrections to the 
Lorentzian line shape. The process dependence of corrections to the Lorentzian line shape in two- 
photon transitions is analyzed, and an order-of-magnitude estimate for off-resonant corrections 
to the peak of the photon scattering cross section is given [see Eq. (|40|) ]. 
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